Unusual thermodynamical and transport signatures of the BCS Bose-Einstein 

crossover scenario below T c 
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In this paper we present predictions for thermodynamic and transport properties of a BCS Bose- 
Einstein crossover theory, below T c , which theory satisfies the reasonable constraints that it yield (i) 
the Leggett ground state and (ii) BCS theory at weak coupling and all temperatures T. The nature 
of the strong coupling limit is inferred, along with the behavior of the Knight shift, superfluid density, 
and specific heat. Comparisons with existing data on short coherence length superconductors, such 
as organic and high T c systems, are presented, which provide some support for the present picture. 
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Within the pseudogap regime of the cuprates, it has 
been widely argued that the excitations of the super- 
conducting state are either fermionic Jl| or bosonic Q in 
character. In this paper we discuss a third scenario (asso- 
ciated with the BCS Bose-Einstein condensation (BEC) 
crossover approach), in which the excitations contain a 
mix of bosonic and fermionic properties . 

The BCS-BEC crossover scheme has been viewed as 
relevant to the cuprates and other "exotic" super- 
conductors where their short coherence lengths, £, nat- 
urally lead to a breakdown of strict BCS theory. The 
BCS-BEC scenario §-|,| owes its origin to Eagles |(J 
and Leggett jll| who proposed a ground state wavefunc- 
tion of the BCS form <3>o = IIk(uk + Wkcj i |cL k ^)|0), which 
describes the continuous evolution between a BCS sys- 
tem, having weak coupling g and large £, towards a BEC 
system with large g and small £. Here itk,fk are the 
standard coherence factors of BCS theory, which are de- 
termined in conjunction with the number constraint. 

The essence of this paper is a characterization of the 
excitations of Wo and their experimental signatures (for 
all T < T c ). New thermodynamical effects stemming 
from bosonic degrees of freedom must necessarily enter, 
as one crosses out of the BCS regime, towards BEC. Here 
we show that the bosonic excitations (which appear in the 
gap equations and which we call "pairons" ) are different 
from the collective phase mode jl^] : they generally have a 
quadratic dispersion similar to that of a quasi-ideal BEC 
system — a consequence of the mean field treatment of 
the pairs which, in turn, is dictated by the general mean 
field character of Wo. These ideas are applied to the 
cuprates and other short £ superconductors. 

We consider fermions, with lattice dispersion ek, (mea- 
sured with respect to the fermionic chemical potential 
/i), and with interaction Vk,k' = SVk^k', where g < 0; 
here tp^ = 1 and (cos k x — cos k y ) for s- and c£-wave pair- 
ing, respectively. We begin by reviewing BCS theory, in 
terms of a somewhat unfamiliar but very useful formal- 
ism which is then generalized |||l4| to the crossover 
problem. For brevity, we use a four-momentum notation: 
K = (k,iw),2 K = T^ ku , etc. We also suppress 9? k 



until the final equations. 

BCS theory involves the pair susceptibility x(Q) = 
J2k G(K)Gq(Q — K), where the Green's function G sat- 



isfies G 1 
tion is 



E, with order parameter A sc and 



A^ c Go(— K). In this notation, the gap equa- 



1+.9X(0) = 0, 



T<T C 



(1) 



At Q = 0, the summand in \ is the Gor'kov U F" 
function (up to a multiplicative factor A sc ) and this 
serves to highlight the central role played in BCS the- 
ory by the more general quantity G(K)Gq{Q — K). Note 
that (for Q ^ 0), xiQ) is distinct from the pair sus- 
ceptibility of the collective phase mode which enters as 
YijAGiK) [G{Q -K) + G(-Q- K)] + 2F(K)F(K - Q)} 
|l2| , [l5| . Here each Gor'kov "F" function introduces one 
GGo, so that the collective mode propagator depends on 
effectively higher order Green's functions than does the 
gap equation. 

The observations in italics were first made in Ref. |i~3| ] 
where it was noted that the BCS gap equation could 
be rederived by truncating the equations of motion so 
that only the one (G) and two particle (T) propaga- 
tors appeared. Here G depends on S which in turn 
depends on T. In general T has two additive con- 
tributions O], from the condensate (sc) and the non- 
condensed (pg) pairs. Similarly the associated self energy 
H E(iT) = £ Q T(Q)G (Q - K) can be decomposed 
into ~E p g(K) + Y, SC (K). The two contributions in E come 
respectively from T SC (Q) = —A^ C 5(Q)/T, and from the 
Q 7^ pairs, with T pg (Q) = g/{l+gx{Q))- In the leading 
order mean field theory E = E sc = — A^ c Go(— K) which, 
from Eq. (1), yields the usual BCS gap equation. 

More generally, at larger g, the above equations hold 
but we now include feedback into Eq. (1) from the finite 
momentum pairs, via T, pg (K) = Y^,q 1~ pg (Q)Go(Q — K) » 
Go(-K)J2 Q T pg (Q) = -A 2 pg G (-K), which defines a 
pseudogap parameter, A pg . This last approximation is 
valid only because (through Eq. (1)) T pg diverges as 
Q — > 0. In this way, T, pg (K) has a BCS-like form, as 
does the total self energy T,(K) = —A 2 Gq{—K), where 



1 



A 2 = A 2 C + A 2 g . Thus, in the present approach, the en- 
ergy gap for single electron excitations reflects the pres- 
ence of both finite center-of-mass momentum pairs as 
well as the condensate. While the structure of the gap 
equation will be seen to be formally identical to that in 
BCS theory, the vanishing of the excitation gap, A, takes 
place at a higher temperature than that at which the or- 
der parameter, A sc , vanishes. The latter defines T c . 



If we now expand 71„ 1 (q, O) w aifi 2 + clqCI 



Bq +iT', we see that the chemical potential of the pairs 
Mpair is proportional to To, and via Eq. (1), precisely zero 
at and below T c . This provides an interpretation along 
the lines of ideal Bose gas condensation. (Here, also, at 
small q, r q — > 0). As g increases, the term ao$l in T~ g Y 
becomes progressively dominant with respect to a\£l 2 . 
For the physically relevant regime of moderate g, we have 
found, after detailed numerical calculations, that a\ may 
be safely neglected. At weak coupling, there is no loss 
of generality in approximating T pg in this more particle- 
hole asymmetric way, since its contribution is negligible. 
In this way, we can write 



^ )fl 1 (q, fi) = a (VL - f2 q + [ipair + ir q 
As a consequence we have 



Q 



(2) 



(3) 



q#0 



We now rewrite Eq. (1), along with the fermion number 
constraint, as 



1 - 2/(£ k ) 2 

2E, ^ k = (K 



E 



i-* + ^/(4) 



(4) 
(5) 



Here f(x) and b(x) are the Fermi and Bose functions and 
Ek = \J(- 2 i + A 2 !/? 2 , is the quasiparticle dispersion. Equa- 
tions (3)-(5) are consistent with BCS theory at small g, 
and with the ground state at all g; in both cases the 
right hand side of Eq. (3) is zero. The simplest physical 
interpretation of the present decoupling scheme is that 
it goes beyond the standard BCS mean field treatment 
of the single particles (which also acquire a self energy 
from the finite q pairs) but it treats the pairs at a self 
consistent, mean field level. 

The dispersion f2 q = q 2 B/ao, as well as the coefficient 
do, are determined by a Taylor expansion of Tsz 1 fl(f | 



T p - X (q,n) 



E 



k 



1 - f(Ek) 



Ek + £k-c 



/(ek-q) 2 

- tlx. 



/(£ k )-/(e k _ q ) 2 



£k-q 



Vk-q/2 



(6) 



For small g and in three dimensions (3D), the poles of T pg 
(at T — 0) occur at f2 = \/3cq, where c is the usual phase 



mode velocity. At moderate g, where the pairons become 
increasingly more relevant, and for quasi-2D dispersion 



£k 



M* /M 



q 2 /2M* 



-q ± /2M]_. Here we find that the ratio 
^ (tj_/t\\) 2 , where t|| and t±_ are the in- and 
out-of-plane hopping integrals, respectively. Numerical 
calculations show that the masses, as well as the residue 
ao, are roughly T independent constants at low T ]l7| ]. In 
the BEC regime at low density and with s-wave pairing in 
a 3D continuous model, M* is 2m e for all T <T C jl6| , as 
found previously 0. The examples in this paper, which 
apply to the fermionic regime, correspond to somewhat 
smaller M«. 

It is important to note that in strictly 2D, the loga- 
rithmic divergence on the right hand side of the pseudo- 
gap equation (3), (which is essentially a boson number 
equation), implies T c = 0, as in an ideal Bose gas. For 
large anisotropy, or small t±, T c cx — l/\a(t±/tu), which 
vanishes logarithmically [H. Finally, since both ao and 
the effective pair mass (tensor) M* are constants at low 
T, Eq. © implies A 2 pg (T) = A 2 (T) - A 2 C (T) cx T 3 / 2 . 
Moreover, because A depends on T only exponentially, 
A 2 C (T) = A 2 (0) - AT 3 / 2 at low T, where A is T inde- 
pendent. 

We now calculate physical quantities such as the mag- 
netic penetration depth (A) and related superfluid den- 
sity (n s ), the Knight shift (K s ) and NMR relaxation 
rate (R s ) using techniques similar to those used to 
study fluctuation effects in normal metal superconduc- 
tors p8| . Here the usual (lowest order) Maki-Thompson 
and Aslamazov-Larkin diagrams are extended to be 
compatible with S and T by applying the generalized 
Ward identity to incorporate the pairon vertex correc- 
tion |l^,[l2|]. Finally, the specific heat C v can be com- 
puted following a similar analysis as for the paramagnon 
problem |]|. 

While in the BCS limit the expressions for A, (or n s ), 
K s and R s contain only the total gap A, here they, in 
principle, depend on both the quasiparticle (via A 2 ) and 

This decom- 



the pairon (via A 2 ) contributions (14 



position leads to a form of "three fluid" model (includ- 
ing the condensate, fermionic quasi-particles and bosonic 
pairons). In the same way, C v can also be decomposed 
into a sum of two contributions corresponding to an ideal 
Bose gas of pairons, with dispersion f2 q , and an ideal 
Fermi gas of quasiparticles, with dispersion E k . The pa- 
iron contributions to n s enter as follows HJl4|]: the gen- 
eral expression is identical to its BCS counterpart, but 
with the overall multiplicative factor of A 2 replaced by 
A 2 C = A 2 — A 2 g . By contrast, for spin-singlet pairing, 
there is no explicit pairon contribution to K s and R s , and 
the corresponding expressions reflect the generalized ex- 
citation gap A, as might have been expected physically. 
The single most important conclusion of this analysis is 
that the presence of low lying pair excitations will intro- 
duce new low temperature power law dependences with 
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FIG. 1. Temperature dependence of (a) Knight shift, K 3 , 
specific heat C v /T (inset), and (b) superfhiid density n s in 
conventional s- and d-wave (BCS) and short £ (PG) .s-wave 
cases, calculated at (n, — g/4£||)= : s-BCS: (0.5, 0.5); s-PG: 
(0.5, 0.7); d-BCS: (0.8, 0.225); d-PG: (0.92, 0.56). 



ideal Bose gas character into physical quantities. Be- 
low, we explore these power laws in the context of highly 
anisotropic 3D, i.e., quasi-2D systems. 

Figures 1(a) and 1(b) present a comparison between 
an s-wave short £ pseudogap (PG) superconductor and 
an s- and d-wave BCS system. It should be noted that 
the short £ superconductors are still far from the BEC 
limit. For the parameters illustrated by the figures, /i 
deviates from Ep by roughly 3%. Here, and through- 
out this paper we take t^/t\\ — 0.01. The main body 
of Fig. 1(a) indicates that the Knight shift (and NMR 
relaxation rate, not shown) at T c are substantially re- 
duced relative to their high T asymptotes, i.e., K n , as 
is illustrated by the solid line (for the PG s-wave). Be- 
cause pairon effects are not explicit, the low T behavior is 
exponentially activated as for the BCS s-wave case, but 
here the ratio A(0)/T c is significantly enhanced over the 
BCS value. Overall, the behavior of R s will yield rather 
similar plots; however, the s-wave BCS limit exhibits the 
well-known Hebel-Slichter peak, which is absent below T c 
for the other two cases. 

In the inset to Fig. 1(a) we plot the behavior of the 
low T specific heat "coefficient", j(T) = C v /T, for the 
same parameters as above. In the short £, quasi-2D case, 
slightly above T = 0, j(T) will appear to be a constant 
7 (T) = 7*, although it vanishes strictly at T = as T 1 / 2 . 
This intrinsic 7* effect, which may have been seen in both 
organic and cuprate (layered) superconductors |^(|, has, 
in the past, been related to extrinsic effects. Figure 1(b) 
plots the normalized superfhiid density n s , or A" 2 , vs 
T/T c , which for the PG (s-wave) case exhibits a T 3 / 2 
dependence. Here A, (unlike C v ), is not particularly sen- 
sitive to the mass anisotropy ratio, and the boson power 
law dependence is more 3D. 

In order to address d-wave effects in short £ supercon- 
ductors, we turn to the cuprates. Note, the dimension- 
less coupling is g/tn. To be consistent with the observed 
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FIG. 2. (a)-(b) Scaling behavior of the T dependence of 
Knight shift, K s , for the cuprates with respect to doping x, 
and (c) doping dependence of "f(T) — 7* + aT and K s at 
T~ . In (a) [and (b)], x varies from 0.05 to 0.2 from top to 
bottom. For comparison, shown in the inset are experimen- 
tal data from Ref. on under- (■), optimally (o, •), and 
overdoped (♦) Bi2Sr2CaCu208+i single crystals. As per the 
experimental convention for <^k, we use 2A(0) in (b). 

metal-insulator transition at half filling (x = 0), we intro- 
duce a hole concentration x dependent renormalization of 
the in-plane hopping integral iy (x) — tgx deriving from 
Coulomb correlations, and presume, in the absence of any 
more detailed information about the pairing mechanism, 
that g is x independent. Our quasi-2D band structure 
is taken from the literature pl[ ; the one free parameter 
—g/4t is chosen (= 0.045) to optimize agreement with 
the energy scales in the cuprate phase diagram. This cal- 
culated phase diagram ||, deduced from Eqs. (3)-(5), can 
be shown to yield reasonable agreement with experimen- 
tal data. Two important points should be stressed: (i) 
the chemical potential ^jl/Ef differs from unity by at most 
a few percent over the entire range of x; (ii) While the 
band mass increases with underdoping due to Coulomb 
effects, the thermodynamically measured mass, obtained 
from, say, the Knight shift K S (T C ) (or specific heat C v /T 
at T~) decreases with underdoping, as a consequence of 
the opening of the pseudogap. [See Fig. 2(c)]. 

Figures 2(a) and 2(b) illustrate the predicted behavior 
of the Knight shift for the cuprates as a function of x and 
T . Because it depends only on (d-wave nodal) quasiparti- 
cle excitations, K s exhibits a scaling with T/A(0) which 
is illustrated in Fig. 2(b) via plots of K s (normalized to 
its high T asymptote K n ) for the entire range of x, and 
for temperatures below each respective T c (x). An alter- 
nate scaling form is shown in Fig. 2(a) where we plot K s 
normalized at T c as a function of T/T c for various x (with 
x increasing from top to bottom). The near-collapse of 
the different x-dependent curves is similar to that found 
in the experimental data p2| ] shown in the inset. The 
normalization factor K S (T C ) for this figure, (which varies 
as the band mass multiplied by T c /A(0)), is plotted as a 
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FIG. 3. Comparison of penetration depth data ^3), AA, 
along a-axis, in nominally pure YBCO6.95 single crystal, with 
different theoretical fits corresponding to BCS d-wave (dashed 
curve) and to BCS-BEC (solid curve) predictions. The cor- 
responding derivatives are plotted in the lower inset. In the 
upper inset are experimental data (AA vs T) for the organic 
superconductor BEDT from Ref. H. 

function of x in Fig. 2(c). Also plotted here are our spe- 
cific heat (C v /T) predictions for the pairon contribution 
to 7* as compared with the usual d-wave quasiparticle 
term |20| ] aT c as a function of x. The pairon term be- 
comes increasingly more important with underdoping. 

In Fig. 3 we present a-axis penetration depth data, 
AA(T), in nominally clean optimally doped YBCO single 
crystal, from Ref. pa], along with d-wave fits to our BCS- 
BEC theory and to a straight line associated with a BCS 
superconductor. Because these two fitted curves are es- 
sentially indistinguishable, in the lower inset we plot the 
slopes dX/dT where the difference between the two sets 
of curves is more apparent. Here it is shown that the low 
temperature downturn of the derivative, seen to a greater 
or lesser extent in all AA(T) measurements, fits our pre- 
dicted T 1 / 2 + const, dependence rather well. This down- 
turn has been frequently associated with impurity effects, 
which yield a linear in T slope for AA at very low T, and 
in this case, provide a poorer fit. While these cuprate 
experiments were performed on a nearly optimal sample, 
the same analysis of an underdoped material yielded sim- 
ilarly good agreement, but with a T 3//2 coefficient about a 
factor of two larger. Future more precise and systematic 
low T experiments on additional underdoped samples are 
needed. Plotted in the upper inset are data |24|] on the 
organic superconductor K-(ET) 2 Cu[N(CN) 2 ]Br (BEDT, 
T c w 11K) which fit a pure T 3 / 2 power law over a wide 
temperature regime; in contrast to the cuprates, there is 
no leading order linear term. At present, there seems to 
be no other explanation (besides the pairon mechanism 
presented here) for this unusual power law at the lowest 
temperatures. 

In summary, within a BCS-BEC crossover theory, 



(based on the Leggett ground state), we find that new 
low T power laws associated with a quasi-ideal gas of 
bosonic pair excitations appear in the thermodynamic 
and transport properties, which may be generally rele- 
vant to short £ superconductors. 
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